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Abstract 

The hadronic charmless B decays into a scalar meson and a vector meson are studied within the 
fi'amework of QCD factorization. The main results are: (i) The decay rates for the /o(980)i^*~ and 
fo{980)K modes depend on the fo — <J mixing angle 9. The experimental measurements can be 
accommodated for 9 ~ 20°. (ii) If the ao(980) is a qq bound state, the predicted branching ratios for 
the channels Oq and af^p^ will be very large, of order 30 x 10^^ and 23 x 10^^, respectively. If the 
observation of or the experimental limit on theses two modes is much smaller than the expectation 
of 25 X 10^^, this could hint at a four-quark nature for the ao(980). (iii) For the ao(1450) 
channels, aQ'(1450)/9^ and ao(1450)/9~ are found to have branching ratios of order 16 x 10^^ and 
22 X 10~^, respectively. A measurement of them at the predicted level will favor the qq structure for 
the ao(1450). (iv) Contrary to the naive expectation that T{B~ — > a^p~) ~ \^{B^ — > Qq p~), we 
found that they have comparable rates due to additional contributions to the a^p~ mode from the Oq 
emission, (v) The predicted central value of B{B Kq (1430)i?!>) is larger than experiment, though 
it can be accommodated within theoretical errors. The decays B A'q(1430)/9 are expected to 
have rates substantially larger than that of B ^ A'Q(1430)7r owing to the constructive (destructive) 
interference between the 04 and ag penguin terms in the former (latter). Experimentally, it is 
thus important to check if the B — > K^p modes are enhanced relative to the corresponding KqTt 
channels. 
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I. INTRODUCTION 



Recently we have studied the hadronic charmless B decays into a scalar meson and a pseudoscalar 
meson within the framework of QCD factorization (QCDF) pj]. It is known that the identification 
of scalar mesons is difficult experimentally and the underlying structure of scalar mesons is not 
well established theoretically (for a review, see e.g. |j]). The experimental measurements of 

B SP will provide valuable information on the nature of the even-parity mesons. For example, 
it was pointed out in [l'] that the predicted fi^ — > a^(980)7r^ and Uq {980)K~ rates exceed the 
current experimental limits, favoring a four-quark nature for the ao(980). The decay B^ — > k~^K~ 
also provides a nice ground for testing the 4-quark and 2-quark structure of the k (or Kq{800)) 
meson. It can proceed through VF-exchange and hence is quite suppressed if the k is made of qq 
quarks, while it receives a tree contribution if the k is predominately a four-quark state. Hence, an 
observation of this channel at the level of > 10~^ may imply a four-quark assignment for the k [i| 

In this work we shall generalize our previous study to the decays B SV {S: scalar meson, V: 
vector meson), motivated by the recent observation of the Kq (1430)</> and fo{980)K*~ modes by 
BaBar: 



B{B- ^ /o(980)K*-;/o(980) ^ vr+vr- 
^ /o(980):^*°; /o(980) ^ ^+^- 
B{B- ^ /o(980)p-;/o(980) ^ ^+^- 
B(^ ^ /o(980)/9°;/o(980) ^ ^+^- 
B(^ ^ /o(980)a;;/o(980) ^ tt+tt- 



= (4.6 ±0.7 ±0.6) X IQ-^ [5], 

= (5.2 ±1.2 ±0.5) X IQ-^ [6], 

= (2.6 ± 0.6 ± 0.9) X 10"^ < 4.3 x 10"*^ [6], 

< 1.9 X 10"^ [7], 

< 0.53 X 10"^ [8], 

< 1.5 X 10"^ [9]. (1.1) 

Recently, the decay B Kq (1430)i?i> has been studied in [IQj] within the framework of generalized 
factorization in which the nonfactorizable effects are described by the parameter A^^^, the effective 
number of colors. The result is sensitive to A'^^^. For example, the branching ratio is predicted to be 
(7.70, 3.95, 1.84) x 10^^ for A'^^^ = (2, 3, 5). Hence, in the absence of information for nonfactorizable 
effects, one cannot make a precise prediction of the branching ratio. A QCDF calculation of this 
and other modes will be presented in this work. 

Since B SP decays have been systematically explored in Q| , it is straightforward to generalize 
the study to the SV modes. In the sector of odd-parity mesons, it is known that the rates of the 
penguin-dominated modes K*tt and Kp are smaller than that of the corresponding Ktt ones by a 
factor of ~ 2. This can be understood as follows. In the factorization approach, the penguin terms 
ag and ag are absent in the decay amplitudes oi B ^ K*Tr, while the effective Wilson coefficients 
a4 and ag contribute destructively to i? — > Kp. In contrast, the tree-dominated pn modes have 
rates larger than that of vrvr with the same charge assignment due mainly to the fact that the p 
meson has a decay constant larger than the pion. We shall see in the present work that the same 
analog is not always true for the scalar meson sector. For example, we will show that the rates for 
Ko°(1430)/9-'° are larger than that of ^o°(1430)7r-'°. 

The layout of the present paper is as follows. In Sec. II we introduce the input quantities 
relevant to the present work, such as the decay constants, form factors and light-cone distribution 
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amplitudes. We then apply QCD factorization in Sec. Ill to study B SV decays. Results and 
discussions are presented in Sec. IV. Sec. V contains our conclusions. The factorizable amplitudes 
of various B — > SV decays are summarized in Appendix A. 

II. INPUT QUANTITIES 

Since most of the essential input quantities are already discussed in [l|, here we shall just 
recapitulate the main inputs. 

A. Decay constants and form factors 

Decay constants of scalar and vector mesons are defined as 

{V{p)\q2lt.qi\0) = fvmve*^, {S{p)\q2l^qi\0) = fsPf., {S\q2qi\0) = msfs- (2.1) 

For vector mesons, there is an additional transverse decay constant defined by 

{V{p,e*)\qcj^,q'\Q) = f^{p^el-p,e;), (2.2) 

which is scale dependent. The neutral scalar mesons a, /o and Oq cannot be produced via the 
vector current owing to charge conjugation invariance or conservation of vector current: 

= //o = /ao = 0. (2.3) 

For other scalar mesons, the vector decay constant fs and the scale-dependent scalar decay constant 
fs are related by equations of motion 

fJ-sfs = fs, with fj.s = — r~\ TT' (^•^) 

where m2 and mi are the running current quark masses and is the scalar meson mass. For the 
neutral scalar mesons /o, and a, fs vanishes, but the quantity fs = fsl^s remains finite. 

In [l| we have applied the QCD sum rule method to estimate various decay constants for scalar 
mesons which are summarized as follows: 

/ao(980)(l GeV) = (365 ± 20) MeV, /,„(980)(2.1 GeV) = (450 ± 25) MeV, 

4(980) (1 GeV) = (370 ± 20) MeV, 4(980) (2-1 GeV) = (460 ± 25) MeV, 

4(1450) (1 GeV) = (460 ± 50) MeV, 4(i450)(2.1 GeV) = (570 ± 60) MeV, 

/i^'{i430)(l GeV) = (445 ± 50) MeV, 4* (1430) (2-1 GeV) = (550 ± 60) MeV. (2.5) 

In [l| we have considered two different scenarios for the scalar mesons above 1 GeV, which will be 
briefly discussed in Sec. IV. The above decay constants for the ao(1450) and Kg (1430) are obtained 
in scenario II. Using the running quark masses given in Eq. ()A13p we obtain the scale-independent 
decay constants: 

/ao(980)± = 1.0 MeV, 4(i450)± = 5.3 MeV, /x*(i430) = 35.9 MeV. (2.6) 
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TABLE I: Form factors for B ^ p, K* , ao(1450), (1430) transitions obtained from the covariant 



light-front model [13 1 



F 


F(0) 


Vymax/ 


a 


b 


F 


F(0) 


V ^ max / 


a 


b 




0.26 


0.68 


1.57 


0.70 


^0 


0.26 


0.35 


0.55 


0.03 


^BK* (1430) 


0.26 


0.70 


1.52 


0.64 


BK* {M30) 
^0 


0.26 


0.33 


0.44 


0.05 




0.27 


0.79 


1.84 


1.28 




0.28 


0.76 


1.73 


1.20 


yBK* 


0.31 


0.96 


1.79 


1.18 




0.31 


0.87 


1.68 


1.08 



For longitudinal and transverse decay constants of the vector mesons, we use (in MeV) 



fp 
Jp 



216 ±3, 
= 165 ± 9, 



J U. 



187 ±5, 
= 151 ±9, 



fK* 

Ik* 



220 ±5, 
: 185 ± 10, 



215 ±5, 
= 186 ± 9 . 



(2.7) 



where the values of fy and fy are taken from 11 1 . 

Form factors for B ^ S,V transitions are defined by [1 

1 



{V{p')\V^\B{p)) 
{V{p')\A^\B{p)) 

{S{p')\A,\B{p)) = 

where = {p + p')^, = 
transition in order for the B 



ruB + my 
i[{mB + my)e*Af^(g 

£* ■ P 



e* -P 



rriB + my 



P„A 



p-^2 



-2my'-fq,[A^^{q^)-A^^{q^, 



Ff'{q') 



+ 



ml 



^%Fi^{q^) 



{21 



(p — p')^. As shown in [l^, a factor of {—i) is needed in B ^ S 
S form factors to be positive. This also can be checked from heavy 



"or B ^ S,V transitions have been evaluated in the relativistic covariant 



quark symmetry 13l |. 

Various form factors 
light-front quark model 

and their momentum dependence is fitted to a 3-parameter form 

.2^ _ ^(0) 



131]. In this model form factors are first calculated in the spacelike region 
ence is fitted to a 3-parameter form 
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a{q^ /my) + b{q'^ /mIY 



The parameters a, b and -F(O) are first determined in the spacelike region. This parametrization is 
then analytically continued to the timelike region to determine the physical form factors at > 0. 
The results relevant for our purposes are summarized in Table HI The form factors for B to /o(980) 



and ao(980) transitions are taken to be 0.25 at g = [l|]. 

We need to pay a special attention to the decay constants and form factors for the /o(980). What 
is the quark structure of the light scalar mesons below or near 1 GeV has been quite controversial. 
In this work we shall consider the conventional qq assignment for the light scalars /o(980) and 
oo(980). In the naive quark model, the fiavor wave functions of the /o(980) and cr(600) read 



a 



-^{uu + dd), 



fo 



ss, 



(2.10) 
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where the ideal mixing for /q and a has been assumed. In this picture, /o(980) is purely an ss 
state. However, there also exist some experimental evidences indicating that /o(980) is not purely 
an ss state (see for details). Therefore, isoscalars cj(600) and /o must have a mixing 

|/o(980)) = \ss)cose + \nn)sme, |o-(600)) = -|ss) sin 6^ + |nn) cos 0, (2.11) 

with nn = {uu + dd) / \/2 . Experimental implications for the fo—cr mixing angle have been discussed 
in detail in [3]. In this work, we shall use 9 = 20°, which is favored by the phenomenological 
analysis of i? — > foK* decays (see below). In the decay amplitudes involving the /o(980) we will 
use the superscripts q = u,d, s to indicate that it is the q quark content of the /o(980) that gets 
involved in the interaction. For example, = ff^ cos 9 and smy/v^. 

B. Distribution amplitudes 

The twist-2 light-cone distribution amplitude (LCDA) $5(2;) and twist-3 ^g{x) and ^5(2;) for 
the scalar meson S respect the normalization conditions 

dx<^six) = fs, C dx^%{x) = [' dx^U^) = fs. (2.12) 
Jo Jo Jo 

In general, the twist-2 light-cone distribution amplitude $5 has the form 

00 



$s(x,/i) = fsQx{l - x) 



m=l 



(2.13) 



where are Gegenbauer moments and Cm^ are the Gegenbauer polynomials. For the neutral 



scalar mesons fo,aQ,a, only odd Gegenbauer polynomials contribute. In p| we have applied the 
QCD sum rules to determine the Gegenbauer moments Bi and B^ (see Table HI]). For twist-3 
LCDAs, we use 

$|(x) = /5, ^Ux) = fs6x{l-x). (2.14) 



TABLE II: Gegenbauer moments Bi and -B3 at the scales /i = 1 GeV and 2.1 GeV (shown in 
parentheses) obtained using the QCD sum rule method [l|. 



State 




Bi 






ao(980) 


-0.93±0.10(- 


-0.64 ± 0.07) 


0.14 ±0.08 


(0.08 ±0.04) 


/o(980) 


-0.78±0.08(- 


-0.54 ± 0.06) 


0.02 ±0.07 


(0.01 ±0.04) 


ao(1450) 


-0.58±0.12(- 


-0.40 ± 0.08) 


-0.49±0.15(- 


-0.29 ± 0.09) 


i^o*(1430) 


-0.57±0.13(- 


-0.39 ± 0.09) 


-0.42±0.22(- 


-0.25 ±0.13) 



For vector mesons, the normalization for the twist-2 function and the twist-3 function is 
given by [3] 

dx^vix) = fv, C dx<^v{x) = 0, (2.15) 
Jo Jo 
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where the definitions for $^(x) can be found in [l5| . The general expressions of these LCDAs read 



= 6x{l - x)fv 



n=l 



and 



2x 



n=l 



(2.16) 



(2.17) 



where Pn{x) are the Legendre polynomials. The Gegenbauer moments a)^ and j_ have been 



studied using t 
at ^ = 1 GeV 



re QCD sum rule method. Here we employ the most recent updated values evaluated 



i6i 



0.03 ±0.02, afj_ = 0.04 ± 0.03, a^ = 0.11 ±0.09 



K* 



a: 



0.15 ±0.07, a: 



0.14 ±0.06, a? = 0.18 ±0.08, 



"2,1 



0.10 ±0.08, 

0.14 ±0.07, (2.18) 



and aY = 0, a]^_|_=0 for V = p,uj, 4>. 



As stressed in [l], it is most suitable to define the LCDAs of scalar mesons including decay 
constants. However, in order to make connections between B SV and B VV amplitudes, 
it is more convenient to factor out the decay constants in the LCDAs and put them back in the 
appropriate places. In the ensuing discussions, we will use the LCDAs with the decay constants 
fsJsJvJv, fp being factored out. 



III. B ^ SV DECAYS IN QCD FACTORIZATION 



We shall use the QCD factorization approach [I5l . Il7l | to study the short-distance contributions 
totheS ^ S'ydecayswithS' = /o (980), ao (980), ao( 1450), i^o*( 1430), and y = p, ET*, cj. In QCD 
factorization, the factorizable amplitudes of above-mentioned decays are summarized in Appendix 
A. The effective parameters af with p = u,c in Eq. ()A8P can be calculated in the QCD factoriza- 
tion approach 17|]. They are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions 
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4{MiM2 



Ci + 



Ni{M2) + 



Ci±l CpOls 



Ci±l 



V,{M2) + —Hi{MiM2)\ + Pf (M2), (3.1 



where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), q are the Wilson coefficients, 
Cf = {N^ - l)/{2Nc) with = 3, M2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi{M2) account for vertex corrections, Hi{MiM2) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 



quark of the B meson and Pi{M2) 
Ni{M2) reads 

N,{M2) 

Note that Ni{M2) vanishes for i = 
equation in Eq. (I2.15p . 



for penguin contractions. The expression of the quantities 
6, 8 and M2 = V, 



0, i 

1, else. 
6,8 and M2 



(3.2) 

V owing to the consequence of the second 
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The vertex and penguin corrections for SV final states have the same expressions as those for 
PP and PV states and can be found in 



15|, [iTI]. Using the general LCDA 



^m{x, jj) = 6x(l — x) 



n=l 



(3.3) 



and applying Eq. (37) in [15] for vertex corrections, we obtain 



nib 



1 



Vi{M) = 12 In — - 18 Sivr + Sivr af 

jJL 2 



11 



20 ^ 



79 _ 2i7r \ 

36 ~ ~3~y 



(3.4) 



for i = 1-4,9,10, 

VAM) = -121n — + 6- - -3i7r 
/i 2 



11 



2>m a 



M 



21 
20 



a 



M 



79 2m 
36 ~ ~Y I 



for i = 5,7 and 



y,(M) 



-6 for M = 5, 

9 - evri for M = V, 



(3.5) 



(3.6) 



for i = 6, 8 in the NDR scheme for 75. The expressions of Vi{M) up to the term are the same 
as that in [l7| . 

As for the hard spectator function H, it reads 



D[MiM2) Jo p Jo ^ Jo 



1 



drj 
fj 



(3.7) 



for i = 1 — 4, 9, 10, where the upper sign is for Mi = V and the lower sign for Mi = S, 



= d7#^ ['-"-^^^ n^MM j 

D[MiM2) Jo p Jo ^ Jo 



1 



"1 dr] 
V 



(3., 



for i = 5,7 and ffj = for i = 6, 8, ^ = 1 — ^ and r] = 1 — ??, $m (^m) is the twist-2 (twist-3) 
light-cone distribution amplitude of the meson M, and 



The ratios r}^ and are defined as 

A, A. 



y _ 2my /^(/x) 
'^>^^^^- rubin) fv ' 



Z)(5y) = ^(0)m|, I)(FS) = ^^(0)^1. 

2m5 fs{p) 



2m| 



r'^(/i) - - 
^ mbip) fs mb{fi){m2{li) - mi{ii)) 



(3.C 



(3.10) 



For the neutral scalars a, fo and Oq, becomes divergent while fs vanishes. In this case one needs 
to express fsfy by fs^y with 



(3.11) 
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Weak annihilation contributions are described by the terms bi, and 6i,EW in Eq. 
have the expressions 



HI) which 



bi 
b2 

&3,EW 
^4,EW 



Cf 
iV2 



ciA\, 



bs 



Cf 
Cf 

m 



Cf r 
Cf 



C3A\ + C5{Al + aI) + NcCqA: 



1/ 



C4^i + C6^^ 



cwA{ + cg^l > 



(3.12) 



where the subscripts 1,2,3 of A^-^ denote the annihilation amplitudes induced from {V — A)(V — A), 
{V — A){V + A) and {S — P){S + P) operators, respectively, and the superscripts i and / refer to 
gluon emission from the initial and final-state quarks, respectively. Their explicit expressions can 
be obtained from A^J" (VV) for the VV case [l8| with the replacements specified in Eq. (|Aip : 



A\ 



A\ 



A\ 



4 



^v{x)^s{y) 
^v{x)^s{y) 



x(l-xy) 


+ 




1 


+ 


1 


x{l-xy) 


xy'' 


1 


+ 


1 


y(l-xy) 


x'-^y 


1 


+ 


1 



x^y 



+ r^r^^vixWsiy) ^ 



for M1M2 
for M1M2 

for M1M2 



VS, 

sv, 
vs. 



2x 



r^^v{x)^%{y) 



xy{l—xy) 
2?/ 



-r^<D|(x)$y (y) ^^^1^ + r}^<^s{x)'^v{y) ^^^^r^ j , 



2(l+p , „S. 
xtp 

2(l+g) 
xip 



r^^>,(x)cl>5(y) ^ + r^^v{x)niy) x^y 



2{l+x) 
2{l+x) 



for M1M2 = 

for M1M2 = VS, 
for M1M2 = 5y, 

for M1M2 = VS, 
for M1M2 = SV, 



-r^<^i.{x)'^y{y)^il^-r^<^six)<^>M 
A{ = a{ = 0, (3.13) 

where f- • • = vros jj^ dxdy, x = \ — x and y = 1 — y. Note that we have adopted the same convention 
as in [ist] that Mi contains an antiquark from the weak vertex with longitudinal fraction y, while 
M2 contains a quark from the weak vertex with momentum fraction x. 

Using the asymptotic distribution amplitudes for vector mesons and keeping the LCDA of the 
scalar meson to the third Gegenbaur polynomial in Eq. (j2.13p . the annihilation contributions can 
be simplified to 

r r / OQ in 



A\{VS) 
AliVS) 

AliVS) 



Bi{3Xa + 4 - tt^) + I^IOXa + 
Bi{Xa + 29 - 3tt^) + B3 (xa + 



18 
2956 



-IT 



100 



-vr 



r^r^XA{XA-2) 
-4r]^XA{XA-2)^ 



V 



9Bi [^XX - 4Xa - 4 + vr^ ) + IOS3 

2\ ~1 



3X1 



19X^ + ^ + 37r2 

D 



Xa - 2Xa + 



IT 



AiiVS) « evra J -3r^^5(^A - 2) 



Bi{6Xa - 11) + ^3 (2OXA - ^) + 4^Ai2XA - 1)| 



(3.14) 
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for Ml Ms = VS, and 

A{{SV) = -AiiVS), AiiSV) = -A\{VS), 

AliSV) = AiiVS), 4{SV) = -AiiVS), (3.15) 

for M1M2 = SV, where the endpoint divergence is defined in Eq. ()3.16p below. As noticed 
in passing, for neutral scalars a, /o and Gq, one needs to express fs'^x •^^^x ™^ fsf^s by fs- 
Numerically, the dominant annihilation contribution arises from the factorizable penguin-induced 
annihilation characterized by ^3. Physically, this is because the penguin- induced annihilation 
contribution is not subject to helicity suppression. 

Although the parameters Oj(i 7^ 6, 8) and a^^sr^ are formally renormalization scale and 75 scheme 
independent, in practice there exists some residual scale dependence in aj(/i) to finite order. To be 
specific, we shall evaluate the vertex corrections to the decay amplitude at the scale // = In 



contrast, as stressed in 17|, the hard spectator and annihilation contributions should be evaluated 



at the hard-collinear scale Hh = y/Jj^^ with ~ 500 MeV. There is one more serious complication 
about these contributions; that is, while QCD factorization predictions are model independent 
in the nii, 00 limit, power corrections always involve troublesome endpoint divergences. For 
example, the annihilation amplitude has endpoint divergences even at twist-2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and posses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow [l^ to parameterize the endpoint divergence Xa = Jq dx/x in the annihilation 
diagram as 

XA = \n(^yi + PAe"t'^), (3.16) 

with the unknown real parameters pA and (pA- Likewise, the endpoint divergence Xh in the hard 
spectator contributions can be parameterized in a similar manner. 



IV. RESULTS AND DISCUSSIONS 

While it is widely believed that the /o(980) and the ao(980) are predominately four-quark states, 
in practice it is difficult to make quantitative predictions on hadronic B SV decays based on 
the four-quark picture for light scalar mesons as it involves not only the unknown form factors and 
decay constants that are beyond the conventional quark model but also additional nonfactorizable 
contributions that are difficult to estimate. Hence, we shall assume the qq scenario for the /o(980) 
and the ao(980). 

For ao(1450)y and i^^^(1430) channels, we have explored in [l| two possible scenarios for 
the scalar mesons above 1 GeV in the QCD sum rule method: (i) In scenario 1, we treat 
K, ao(980), /o(980) as the lowest lying states, and Kq{14:30), oq (1450), /q (1500) as the corresponding 
first excited states, respectively, where we have assumed that /o(980) and /o(1500) are dominated 
by the ss component and (ii) we assume in scenario 2 that Kg (1430), oq (1450), /o (1500) are the 
lowest lying resonances and the corresponding first excited states lie between (2.0 ~ 2.3) GeV. 
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TABLE III: Branching ratios (in units of 10^^) of B decays to final states containing a scalar 
meson and a vector meson. The theoretical errors correspond to the uncertainties due to (i) the 
Gegenbauer moments -61,3, the scalar meson decay constants, (ii) the heavy-to-light form factors 
and the strange quark mass, and (iii) the power corrections due to weak annihilation and hard 
spectator interactions, respectively. The predicted branching ratios oi B ^ /o(980)K*, /o(980)7r 
are for the /o — o" mixing angle 9 = 20°. For light scalar mesons /o(980) and ao(980) we have 
assumed the qq content for them. The scalar mesons ao(1450) and Kq{14:50) are treated as the 
lowest lying scalar states, corresponding to scenario 2 as explained in Appendices B and C of jlj. 
Experimental results are taken from Eq. (jl.ip . We have assumed i3(/o(980) — > vr^vr^) = 0.50 to 
obtain the experimental branching ratios for /o(980)y. 



Mode 



Theory Expt 



Mode 



Theory Expt 



/o(980)K*' 
/o(980)pO 
/o(980)a; 
a^i980)K*- 
a[J(980)K*° 
a^{980)p- 
Oq (980)p+ 
a[}(980)p° 
ag(980)a; 
a^(1450)is:*- 

afj'(1450)p- 
Oq (1450)p+ 
a[](1450);o0 
a[](1450)w 

i^o"(1430)p+ 
K*o\U30)p^ 36.8 
:^o°(1430)a; 



B- 
B 

B- 
B- 
B- 
B^ 
B 

B~ 
B- 
B- 
B- 
B- 

B- 

B- 
B~ 
B- 



f 0(980) K*- 
fo{980)p- 

a^{980)K*~ 
ao{980)K*° 
aS](980)p~ 
Oq (980)p° 
(980)cj 



aS](1450)K*- 
ao(1450)K*° 
a[](1450)p- 
(1450)p° 
Oq (1450)u; 

K*-{U30)cl) 

Ao*-(1430)p° 
i^o"(1430)u; 



7 /i+o. 4+0. 2+7.2 

' -^-0.4-0.2-2.9 
1 q+0. 1+0.4+0.1 
-■-•"J-0. 1-0.3-0.1 



n q+0.1+0.4+5.6 
^■"-•-0.1-0.4-1.3 
1+0.2+1.1+4.7 
0-J^-0. 2-1.0-2.1 
+0.6+3.6+3.9 
0.5-3.3-2.6 
+0.8+0.1+2.3 
-"-0.6-0.1-1.0 
1 n+0-3+0.0+0.3 
-'-•'-'-0.2-0.0-0.3 



23.0 
3.0 



-0.2- 
-0.2- 



-0.44 
-0.3- 



1.7 



^'-'- '-4.6- 



66.2 
21.0 



10.4 ±2.6 
< 3.8 



12.1 
1.3 



2.2 



7 0+0.5+1. 2+23. 6 
'•°-0.5-l.l- 4.6 
91 7+0.2+3.4+4.2 
' -0.2-3.1-3.9 
rr q+1. 6+0.1+7.6 
"J-"J-1. 3-0.1-2. 4 
+0.6+0.0+1.3 
-0.5-0.0-1.0 



1.6+52.1 
1.6-10.1 
25.0+2.8+70.8 
19.5-2.4-26.3 
+7.3+1.2+29.4 
5.9-1.1-10.1 
-iR 1+4.9+0.7+22.5 
-'-'-'•-'--4.0-0.6- 8.3 



f- ^+0.4+0.3+7.0 
'-'•^-0.4-0.3-2.6 
r, ni +0.00+0.00+0.02 
^•^-'--0.00-0.00-0.01 
hO. 02+0. 00+0.02 
-0.01-0.00-0.02 
-0.1+0.7+3.5 
-0.1-0.6-1.4 
-+0.2+0.5+6.5 
-0.2-0.4-1.8 
-2.4+5.9+6.0 
-2.1-5.3-4.5 
hO.02+0.00+0.30 
'-0. 02-0. 00-0. 07 
2 g+0. 3+0.0+0.2 



0.06] 

3.9^ 
2.5^ 
29.5^ 
0.09^ 



1.3 



-0.3- 

-0.2H 
-0.2- 



0.0- 
0.0- 
0.0- 



4.7 



16.0 



16.4 



5.2 ±2.2 < 8.6 

< 1.06 

< 3.0 



0.2 
-0.1 
■0.1 



+0.3+0.7+19.1 
-0.3-0.7- 2.3 
r, t;-l-0.4+0.5+13.6 
^•O—n -^—n A— 2.2 



0.3-0.4- 
+2.8+4.3+20.8 
2.4-3.8-10.4 
9 9+1.0+0.0+1.3 
^•^-0.7-0.0-0.5 
q 1+0.7+0.0+1.9 
"-'•-'--0.6-0.1-1.2 
1 0+0.5+0.0+1.3 
-•-•0-0.4-0.0-0.8 
+6.1+1.6+51.6 
4.6-1.5-10.1 



cri r,-l-16.1-l-l-4+68.6 
^-'-•'-'-13.1-1.2-23.8 
+14.3+0.9+23.4 
11.0-0.7- 9.1 
1 c f;+4.4+1.0+14.6 
-'-"-'•'J-3.7-0.8- 5.3 



4.6 ±0.9 



Scenario 2 corresponds to the case that light scalar mesons are four-quark bound states, while all 
scalar mesons are made of two quarks in scenario 1. We found that the predicted ao(980)i^ and 
aQ'(980)7r~ rates in scenario 1 will be too large compared to the current limits if the ao(980) is 
a bound state of the conventional qq quarks. This means that the scenario 2 in which the scalar 
mesons above 1 GeV are lowest lying qq scalar state and the light scalar mesons are four-quark 
states is preferable. Therefore, we shall use scenario 2 when discussing ao(1450) and iiro(1430) 
mesons. 

The calculated results for the branching ratios of ^ SV are shown in Table IIIIl In the 
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table we have included theoretical errors arising from the uncertainties in the Gegenbauer moments 
Bi^3 (cf. Table [n|) . the scalar meson decay constant fs or fs (see Eq. (|2.5|) ). the form factors 
pBP,BS ^ the quark masses and the power corrections from weak annihilation and hard spectator 
interactions characterized by the parameters Xa and Xh, respectively. For form factors we assign 
their uncertainties to be (5F^^'^^(0) = ±0.03, for example, ^^^^(O) = 0.35 ± 0.03 and Fq^*°{Q) = 
0.26 ± 0.03. The strange quark mass is taken to be ms(2.1 GeV) = 90 ± 20 MeV. For the quantities 
Xa and Xh we adopt the form (|3.16|) with pa,h ^ 0.5 and arbitrary strong phases 4>a,h- Note 
that the central values (or "default" results) correspond to pa,h = and (pA,H = 0. 



1. B ^ fo{980)K* and oo(980, 1450)i^* decays 

The penguin-dominated B /o(980)i^* decay receives three distinct types of factorizable 
contributions: one from the K* emission, one from the /o emission with the ss content, and the 
other from the /o emission with the nn component. ^ In the expression of S — > foK* decay 
amplitudes given in Eq. ()A8p . the superscript u of the form factor Fq ° reminds us that it is the 
u quark component of /o involved in the form factor transition. In contrast, the superscript q of 
the decay constant fj^ indicates that it is the qq quark content of /o responsible for the penguin 
contribution under consideration. Note that except for the penguin operators Oq and Os, the /o 
emission amplitudes induced from other four-quark operators contain a vanishing /q decay constant. 
However, it is compensated by the ns term in the twist-2 LCDA of the scalar meson so that the 
combination fj^ A*/o = ffo becomes finite. 

In the extreme case that the /o(980) is made of ss quarks or nn quarks, the branching ratio of 
fo{980)K*- is given by 

BiB ^/o(980)K for/o(980)=^n. ^'''^ 

In general, B{B — > /o(980)i^*) depends on the mixing angle 9 of strange and nonstrange components 
of the /o(980). We found that the experimental data can be accommodated with 6 being in the 
vicinity of 20° (see Fig. [1]). The charged and neutral modes of /o(980)-fir* are expected to have 
similar rates, while experimentally their central values differ by a factor of 2. This discrepancy 
needs to be clarified by the future improved measurements. 

In order to compare theory with experiment for B — > fo(980)K* , we need an input for 



;S(/o(980) —>■ vr+TT ). To do this, we shall use the BES measurement [20|] 



r(/o(980) ^ Ti) + r(/„(980) ^ KK) 

Assuming that the dominance of the /o(980) width by vrvr and KK and applying isospin relation, 
we obtain 

S(/o(980) ^ TT+TT-) = 0.50;°:°^ , e(/o(980) ^ K+K') = 0.1251°:°^^ . (4.3) 



^ In our previous work for B SP decays [l|, we did not take into account the contributions from the /o 
or the neutral ao emission induced from the four-quark operators other than Og and Og (see also IS]). 
Corrections will be published elsewhere. 
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0° 90° 180° 0° 90° 180° 

e e 

FIG. 1: Branching ratios of B' fo{980)K*- and B° /o(980)i^*° versus the mixing angle 9 
of strange and nonstrange components of /o(980). For simphcity, theoretical errors are not taken 
into account. The horizontal band shows the experimentally allowed region with one sigma error. 

Hence, we use ;S(/o(980) tt+tt") = 0.50 to determine the absolute branching ratio for B 
M980)K*. 

For aoK* decays, they have similar branching ratios as the corresponding aoK channels 



2. B ^ /o(980)p and ao(980, 1450)p decays 

The tree dominated decays B ao(980)/3, /o(980)p are governed by the B ^ and B ^ Jq 
transition form factors, respectively. The fop rate is rather small because of the small uu component 
in the /o(980) and the destructive interference between 04 and ag penguin terms. The fop^ and 
fou> modes are suppressed relative to fop^ by at least a factor of ^|a2/aip. 

The decay B^ (^oP~ has a rate much larger than the Oq p"*" one because the factorizable 
amplitude of the former (latter) is proportional to fp {fa) and the decay constant of the charged 
oo is very small. We also notice that the predicted aop~ rates are much larger than that of oovr^ 
for two reasons. First of all, the p meson decay constant is bigger than that of the pion, fp ^ f^^. 
Second, the destructive interference between the 04 and ag penguin terms is less severe for qqp as 

~ 2.4r^. Contrary to the naive anticipation that T{B~ AqP") ~ ^T(B^ a^p~), we found 
that they have comparable rates due to additional contributions to the agp" mode from the ag 
emission. 

In the sector of the aQ(980)7r channels, we have argued before that the fact that the experimental 
limits for the apVT and QqK modes are smaller than the theoretical expectations favors a four-quark 
nature for the ao(980) [1]. Here we also suggest that if the observation of or the experimental limit 
on the decay mode a|j'(980)p~ (ao(980)/9~) is much smaller than the expectation of 25 x 10^^, this 
could indicate a four-quark structure for the ao(980). 

Recently, the isovector scalar meson ao(1450) has been confirmed to be a conventional qq meson 
in lattice calculations 



22I. l23l. \24. l25l. l26l]. Hence, the calculations for the ao(1450) channels should 
be more trustworthy. Our results indicate that a,|(1450)p^ and ag(1450)/j~ have large branching 
ratios, of order 16 x 10^^ and 22 x 10^^, respectively. A measurement of them at the predicted 
level will reinforce the qq nature for the ao(1450). 
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3. B ^ i^^(143O)0 and K^{im)p decays 

For i^fj (1430)0 channels, the central value of the predicted ^o(1430)(^) = 

X 10"'^ is larger than the experimental value of (4.6 ± 0.9) x lO"'', though they 
are consistent within theoretical uncertainties. This mode was measured by BaBar [5|] using the 



LASS parametrization to describe the (Ktt)'q' amplitude. However, as commented in 27|], while 
this approach is experimentally motivated, the use of the LASS parametrization is limited to the 
elastic region of M{Ktt) < 2.0 GeV, and an additional amplitude is still required for a satisfactory 
description of the data. Therefore, it will be interesting to see the Belle measurement for Kq{1430)4> 
modes. 

Theoretically, the Kq{1430)p rates are expected to be substantially larger than that of the 
Kq{1^30)tt ones since the penguins terms 04 and ag contribute constructively to the former and 
destructively to the latter. However, as shown in , our predicted central values for the branching 
ratios of K'^I'tt^ and Kq^tt^ are too small by a factor 3 ~ 4 compared to experiment. ^ It appears 
that one needs sizable weak annihilation in order to accommodate the KqTt data. In this work, we 
found large rates for K*q p~~'^ and K^^p^ even in the absence of weak annihilation contributions. 
Experimentally, it should be relatively easy to search for those -fCo(1430)yO modes to see if they are 
enhanced relative to their counterparts in the K^n sector. The branching ratios for the Kq{\4:'3^)u 
modes are predicted to be of order 1.5 x 10~^. 



V. CONCLUSIONS 



We have studied the hadronic charmless B decays into a scalar meson and a vector meson within 
the framework of QCD factorization. The main results are: 

• The decay rates for the /o(980)-R'*~ and /o(980)i^ modes depend on the mixing angle 6 of 
strange and nonstrange components of the /o(980). The experimental measurements can be 
accommodated for 9 ~ 20°. 

• If the ao(980) is a qq bound state, the predicted branching ratios for the channels a^/9~ and 
OqP" will be very large, of order 30 x 10~^ and 23 x 10~^, respectively. If the observation 
of or the experimental limit on theses two modes is much smaller than the expectation of 
~ 25 X 10~^, this could hint at a four-quark nature for the ao(980). 



^ Recently, the authors of [28| claimed that the decay rates for the Kq-k" and Kq~ 
accommodated in the pQCD approach. It is not clear to us what is the underlying reason for the dis- 
crepancy between our work and 



. However, we have just performed a systematical study of charmless 
3-body B decays based on a simply generalized factorization approach ^2!|] . We consider the weak process 
B iirQ(1430)7r followed by the strong decay Kq Kn and reach the same conclusion as 0, namely, 



-rT*0 



the predicted tt and 7r+ rates are too small compared to the data. 
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• For the ao(1450) channels, aQ"(1450)/9~ and ao(1450)/3~ are found to have branching ratios 
of order 16 x 10^^ and 22 x 10~^, respectively. An observation of them at the predicted level 
will favor the qq structure for the ao(1450). 

• Contrary to the naive expectation that r{B~ — aQp") ~ ^^{B^ — 0'oP~), we found that 
they have comparable rates due to additional contributions to the aQp~ mode from the Qq 
emission. 

• The predicted central value of B{B — > Kq (1430)(^) is somewhat larger than experiment, 
though it can be accommodated within theoretical errors. The decays B Kq(1430)/9 
are expected to have rates substantially larger than that of B ^ iCQ(1430)7r owing to the 

constructive (destructive) interference between the 04 and penguin terms in the former 
(latter). Experimentally, it is thus important to check if the B K^p modes are enhanced 
relative to their counterparts in the K^t: sector. The branching ratios for the -ftro(1430)a; 
modes are predicted to be of order 1.5 x 10~^. 
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APPENDIX A 



The B — > SV (VS) decay amplitudes can be either evaluated directly or obtained readily from 
B — > VV amplitudes with the replacements: 

$y(:c)^%(x), ci>^(x) ^ f fv^fs, fv^-fs, ^x^-^x- (Al) 

As stressed in the main text, we use the LCDAs with the decay constants being factored out. Since 
the VV channels have been studied in details in 18|], we may use them to obtain the B — > SV 



amplitudes. In 18|], the factorizable longitudinal B VV amplitude reads (apart from the effective 
Wilson coefficients) 

Av,V2 = i^2fv,A^^'iml^)mBPc, (A2) 

where use has been made of the replacement mye* -pB ^ tt^bPc with pc being the cm. momentum. 
Since the definitions for the decay constant fy and the form factor in 

{V\V^\{)) = -ifvmve% {V{p')\A^\B{p)) = 2mv ^ q^.A^'' {q^) + • • • (A3) 

are different from ours [see Eqs. (j2.ip and (j2.8p ]. the replacements ()Aip need to be modified 
accordingly. The B VS amplitude is obtained from the replacements: 

For B SV amplitudes, the replacements are 

fv.^ifs, fv.^ifv, A^^i^-iFf^, r^^^r^. (A5) 

From (IA4j) and (jASp we obtain the factorizable B SV and VS amplitudes 

^ .GFj2fvF{'^{ml)mBPc; for M1M2 = 5^, 
^M,M, '^\_2fsA^^iml)mBPc; for M1M2 = ^5. ^ ^ 

The coefficients of the flavor operators for 5*^ can be obtained from the VV case [l8^ and 
they read 

ai{Mi,M2) = ai{Mi,M2), a2{MiM2) = a2{MiM2) , 
al{MiM2) = af(MiM2) + af(MiM2), 
al^^{MiM2) = ag(MiM2) + a?(MiM2), 

r al{MiM2) - rlal{MiM2); for M1M2 = SV , 



al{MiM2) 



(M1M2) + r^af (M1M2); for M1M2 = VS, 

a?o(MiM2) - r^4(MiM2); for M1M2 = 5^, 
a?o(MiM2) + r|ag(MiM2); for M1M2 = VS. 



<Ew(^i^2) = <! .P^(M,M2) +r|ag(- - - - - - ^^^^ 

Applying the replacement ()Aip and Eq. ()A7p to the -B — > amplitudes 



m 



181], we obtain the following the factorizable amplitudes of the decays B 
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{fo,aQ)K*, Jq{p,lo), aQ{p,uj), aoK*, K^{(f), p,uj): 
.Gf 



r,=U.C [ 



p=u,c 

- [ a2SP + 2(a3 + as) + ^(07 + ag)) 2/% A^""' {m%)mBPc 



ff^{b2SP + 63 + h,Ew) fuK* + ff^ {^2^1 + 63 + 63,Ew)i^* fs 



V 2 p=„^c I V Z Z ) jdf^. 



- [a3 + al + a5 + {al- -ag)f^° - (a7 + ag + a?o) ) 2f% A^^' (m^Jmep, 



>^ 2 



a25P + 2(a3 + ag) + ^(07 + ag) ) 2/% A^""* {m%)mBP 



fsfK 

.Gf 



p=u,c L "° 
- (a25S)x*ao Vao^O^' {rnl^)rnBPc - fBfK*fao{b2SP + h+ ^3,Ew)„oK* I' 

V 2 p=„^c IV Z Z / a^K* 

X 2/k* «° {mj,. )mBPc - Ib/k* fao (62^^ + ^'3 + h,Ew) aoK* | > 

X 2fK*Ff'^{mj^.)mBPc - fBfK^faoih - ^b3,Ew)aoK- 

A(b' ap<*') = ^^ E - (^^4 - r^K ' Ul) ' Uo) 2fK*Ff^%mh)mBPc 

^ p=u,c IV Z Z JaoK* 



- (a2(5S)if*ao 2/ao-4o {m^JmBPc - fBfK*fao{ - h + 2^3,Ew)aoK* |, 
p=u,c I •'0'' 

X 2/pFf ^° {ml)mBPc - {(^26^ + 2(0,3 + ag) + 04 + (a^ - ^al)f{° 
+ ^(07 + ag - aio))^^^^2//^ AQ^{m%)mBPc 
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- fBfpfl ^{hSP + 63 + &3,EW) fup + {h^P + h + 63,Ew)p/. 

AiB" /oP°) = ^^Y: a(^) ( (a,6^ - < + r^(a^ - U) + ^(af + a?) + ^a?, ^ 

X 2fpF^^°{rnl)mBPc + (as^S + 2(a3 + 05) + 04 + (a| - ^af)f^° 

+ ^(a? + 09 - aio)) ^,„2/]^ AQ''{m%)mBP. 

— fsfpffQ 

p=u,c L /o 

X 2f^Ff^^ml)mBPc - {0,26^ + 2(03 + 05) + 04 + (af - ^af)f^° 



pfS 

- - 1- 3- - - 1- 3- 

{hSP -h + -63,EW + 2^4,EW) fdp + {bi5P - 63 + 2^3,EW + -b4,Ew)pfd 



(5l<5S + ^3 + 264 - 2^3,EW + 2^4,EW) /d^ 



OOP 



+ (^1<^S + ^3 + 264 - ^63,EW + ^^'4,Ew)a,/d 



= E ^i'^ I ("i-^S + <- r'M + «I) + <o 

V ^ p=u,c L 

X 2/pFf {ml)mBPc - fsfpfao \{h + 64 - ;^&3,EW - ;^kEw! 



+ (6l<5S + &4 + 64,Ew)p„o 
V p=u,c K 

11 

X '^faoAo''{mp)mBPc - fsfpfao {h + b4- -63,EW - 2^4,EW; 



aop 



pao 



aop 
ao 



+ {biSP + b4 + b4,Ew)aopi 

A{B-^a^,p-) = E |(ai<5S + a^-r^K + af) + a?o 

p=u,c (. 

X 2fpF^-'{ml)mBPc - {a^Sl - 04 - (a| - ^af)f° 

3 1 
+ 0(^7 + 09) + T^aio) 2faoA^''{ml)mBPc 

Z L ' pao 

- fsfpfao [{hSP + 63 + fe3,Ew)„op - (&25^ + &3 + h,Ew)pao] | ' 

^(S-^ao-p°) = E |-(ai<^S + aI + r^°K + aD + afo) ^ 2/„„Ao^''(m2jmBPc 

p=u,c I. '''''' 
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+ 



113 

a2SZ -al + rP{al - -af) + -afo + -{ag + 07) 



OOP 



- fsfpfao [{b2SP + 63 + b3,Ew)pao ' (^2<5^ + &3 + &3,Ew)„op 



p 



aop 



3, _ , 1 



1 



+ /s/p/ao [(&1<^J^ + ^3 + 264 - ^(^a.EW - b4,Ew))aop 
1 



+ (&1<^^ + ^3 + 264 - 2(^3,EW - fe4,Ew))paoJ 



+ 



a2(5g + 2(03 + 05) + al- r'^ial - -af) 



- ^afo + ^(ag + 07)] 2f^Ff"-%mi)niBPc 

- /B/a;/ao [{b2Sf, + 63 + &3,Ew)a;ao + (&2(5^ + bs + b3,Ew)aou, 

A(^^a',u) = -i^ E Af|(a2<5S + 2(a3 + a5) + a^-r-K-^ai) 

- Jaio + J(a9 + a7)) 2f^Ff''°{ml)mBPc 

- (a2(5S + 2(a3 + a5)+a4 + (ag-^af)f^o 



+ ^(07 + 09 - aio)),,__'^faoAo'^{mljmBPc 



UJOO 



- r - - 1- 3- 

- fsfajfao ( - blSP + 63 - -63,EW - 7;b4,Ew)aoiM 



- - 1- 3- 
+ ( - biSf, + 63 - -63,EW - 2^4,Ew; 



A{B- K*-^) = ^% E I («3 + < + a,- - lag) - l(a7 + ag + a?o)) 

X 2f^Ff^° {ml)mBPc - fBUfK*{b2S^ + 63 + &3,ew)k*</,|> 

= E \ f«3 + a^ + a5 - r^ial - Ja^) - ^(07 + ag + a?o)) 

V 2 I V 2 2 y *^ 

< 2/0Ff ^° {ml)mBPc - IbUIk^ (63 - ^&3,Ew)i^*0| , 
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.Gp 



p=u,c y 



K, 



no 



*)mBPc - fBfpfK*{b25P + 63 + 63,Ew)pii-. 

A{B-^K*-p') = ^^ E A(^)(-(aidS + < + r-?(a^ + aD+<o) 

p=u,c ( '^0 



02(5^ + -(ag +07) 



p=u,c { 



.Gf 



p=u,c 



-«4 - (^6 - ^«8) + ^^10 



fsfpfR'^i - ^3 + 2^3,Ew)pii'*| 



^2(5^ + -(09 + ar) 



2fpFf'^\ml)mBPc 



.Gp 



E A?^- 



p=u,c 



a2SP + 2(03 + 05) + -(09 + ay) 



K* 



fBfufK*{b25l + 63 + bs.Ew)^/^. k 



p=u,c 
K, 



a26^ + 2(03 + as) + -(a9 + ay) 



- (a2 + V(a^ 



2^Ff'^«(m2)mBPc 
|«8)-^«io) ^2fK'^AQp{m\,)mBPc 



fBfujfK'ibs - 2^3,Ew)a;i^*|) 



where Ap'^^ = VpbV* with q = d,s and 



2m2± 



J- 2m f „o 2m„o „± 
?t/q^^^ _ JO f = r (;/) = 

^ mb(/i)' ^ r?T.;,(/x)' ^ mb{p){md{fJ,) - rriuip,))' 



(A9) 



In Eq. (|A8p . we encounter terms such as aj/jg, which appears to vanish at first sight as ff^ = 0. 
However, when fj^ combines with /ij^ appearing in the twist-2 LCDA of the scalar meson [see Eq. 
(j2.13p ]. it becomes finite, namely, //oAi/o = //o- Therefore, the effective Wilson coefficients aj in 
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Eq. (jASP are defined as aifig^ and they can be obtained from Eq. (j3.ip by retaining only those 
terms that are proportional to fis- Specifically, 



Ci±l CpOis 



Vi{M2) + —Hi{MiM2)\ + Pf (M2) 



(AlO) 



The LCDA of the neutral scalar meson in the bar quantities, Vi{S), Pi{S) and Hi{Ml, M2) is 
replaced by $5 which has the similar expression as Eq. ()2.13p except that the first constant term 
does not contribute and the term fs ij-s is factored out 



^six,fi) = 6x{l-x) BMC^\2x-l). 



(All) 



m=l 



InEq. (lATnl) . 



^ - Sivr) Bf + (§-^^Bl + ■■■■ for i = 1 - 4,9, 10, 
^-3m)i?f-(i-f^)i?f + •••; fori = 5,7. 



(A12) 



0; 



for i = 6, 8. 



The annihilation terms bi have the same expressions as Eq. ()3.12p with and iJ-sBi replaced by 
and Bi, respectively. 

For the CKM matrix elements, we use the Wolfenstein parameters A = 0.806, A = 0.22717, 
p = 0.195 and f/ = 0.326 2l|. For the running quark masses we shall use 

mbirrib) = 4.2 GeV, mfc(2.1 GeV) = 4.95 GeV, mfe(l GeV) = 6.89 GeV, 

rricimb) = 1.3 GeV, mc{2.1 GeV) = 1.51 GeV, 

m,(2.1GeV) = 90MeV, m,(lGeV) = 119 MeV, 

mrf(lGeV) = 6.3 MeV, m„(lGeV) = 3.5 MeV. (A13) 
The uncertainty of the strange quark mass is specified as 771^(2.1 GeV) = 90 it 20 MeV. 
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